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1 Introduction 

Experimental observation |l[ of the inclusive and exclusive radiative decays B — > Xg'j and 
B -^ K'y stimulated the study of rare B decays on a new footing. These Flavor Changing 
Neutral Current (FCNC) b —^ s transitions in the Standard Model (SM) do not occur at 
the tree level and appear only at the loop level. Therefore the study of these rare B- 
meson decays can provide a means of testing the detailed structure of the SM at the loop 
level. These decays are also very useful for extracting the values of the Cabibbo-Kobayashi- 
Maskawa (CKM) matrix elements 0], as well as for establishing new physics beyond the 
SMi. 

Currently, the main interest on rare i?-meson decays is focused on decays for which 
the SM predicts large branching ratios and can be potentially measurable in the near 
future. The rare B -^ Ki'^i~a.nd B — > K*£^£^ decays are such decays. The experimental 
situation for these decays is very promising [Q, with e"'"e~ and hadron colliders focusing 
only on the observation of exclusive modes with i = e, ^ and r final states, respectively. 
At quark level the process b — >■ si~^i~ takes place via electromagnetic and Z penguin and 
W box diagrams and are described by three independent Wilson coefficients C7, Cg and 
CiQ. Investigations allow us to study different structures, described by the above mentioned 
Wilson coefficients. In the SM, the measurement of the forward-backward asymmetry and 
invariant dilepton mass distribution in 6 — > qC^Cr {q = s, d) provide information on 
the short distance contributions dominated by the top quark loops and are essential in 
separating the short distance FCNC process from the contributing long distance effects |Q 
and also are very sensitive to the contributions from new physics 0. Recently it has been 
emphasized by Hewett [0 that the longitudinal lepton polarization, which is another parity 
violating observable, is also an important asymmetry and that the lepton polarization in 
b -^ s£~^i~ will be measurable with the high statistics available at the B-factories currently 
under construction. However, in calculating the Branching ratios and other observables 
in hadron level, i.e. for B -^ Ki^i^decaj, we have the problem of computing the matrix 
element of the effective Hamiltonian, 'Heff-, between the states B and K. This problem is 
related to the non-perturbative sector of QCD. 

These matrix elements, in the framework of different approaches such as chiral theory 
0], three point QCD sum rules method ||^, relativistic quark model by the light-front 
formalism |l^, |ll|, have been investigated. The aim of this work is the calculation of these 
matrix elements in light cone QCD sum rules method and to study the lepton polarization 
asymmetry for the exclusive B — > fC£"'"£~decays. 

The effective Hamiltonian for the b -^ st^l' decay, including QCD corrections [12-14] 
can be written as 

K// = -^Vt^VL E C^(/^)0.(/x) , (1) 

v2 i=i 

which is evolved from the electroweak scale down to ;U ~ ?7i(, by the renormalization group 
equations. Here Vij represent the relevant CKM matrix elements, and Oj are a complete 
set of renormalized dimension 5 and 6 operators involving light fields which govern the 
6 — i> s transitions and Ci{fi) are the Wilson coefficients for the corresponding operators. 
The explicit forms of Ci(/x) and Oj(/x) can be found in [12-14]. For b — ^ si~^i~ decay, this 
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where p"^ is the invariant dilepton mass, and L{R) = [1 — (+)75] /2 are the projection 
operators. The coefficient Cg (/i, p'^) = Cg{fi) + Y{^, p^), where the function Y contains 
the contributions from the one loop matrix element of the four-quark operators, can be 
found in [12-14]. Note that the function Y{fi, p^) contains both real and imaginary parts 
(the imaginary part arises when the c-quark in the loop is on the mass shell). 

The B — »• ir£+£~ decay also receives large long distance contributions from the cascade 
process B — > KJ/iplip') — > Ki^i~. These contributions are taken into account by intro- 
ducing a Breit-Wigner form of the resonance propagator and this procedure leads to an 
additional contribution to Co of the form lITsi 
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As we noted earlier, in order to calculate the branching ratios for the exclusive B — »• 
i^£+£~decays, the matrix elements (i^|s7^(l — '~^^)q\B) and {K\sia ^yp'' {1 + '-^^)q\B) must 
be calculated. These matrix elements can be parametrized in terms of the formfactors as 
follows (see also [01): 



{K{q)\-s^,{l-^,)q\B{p + q)) = 2q,f^{p')+[f+{p') + r{p') 

IAp') 
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niB + rriK 



(3) 
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where p + q and q are the momentum of B and K and P^ = {p + 2g) . 



2 Sum rules for the B ^ K transition formfactors 

According to the QCD sum rules ideology, in order to calculate these formfactors, we 
start by considering the representation of a suitable correlator function in hadron and 
quark-gluon languages. For this purpose, we consider the following matrix elements of the 
T-product of two currents between the vacuum state and the /T-meson: 



n«(p,g) = I d'x e'^%K{q)\T{s{x)^,{l-^,)h{x)mi^,q{m^) 



nf (p,g) = I d'x e'P^K{q)\T{s{x)za,,p%l + j,)b{x)b{0)t^,q{0)}\0) 



(5) 
(6) 



where q is /r-meson momentum and p is the transfer momentum. 

The hadronic (physical) part of eqs.(5) and (6) is obtained by inserting a complete set 
of states including the 5-meson ground state, and higher states with 5-meson quantum 
number: 
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where P^ = {p + 2q)^. Then, for the invariant amphtudes Fj, one can write a general 
dispersion relation in the B meson momentum squared, {p + qY, as: 



F.{pMp+qr))= r^ 

J mi S — 
where the spectral densities are given by 
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The first terms in eqs.(lO), (11) and (12) represent the ground state 5-meson contribution 
and follow from eqs.(5) and (6) by inserting the matrix elements given in eqs.(3), (4) and 
by replacing 

fsml 



{B\bt^,q\0) 



rrib 
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in which we neglect the mass of the light quarks. In eqs.(10)-(12), p^{p^.,s) represent the 
spectral density of the higher resonances and of the continuum of states. In accordance with 
the QCD sum rules method we invoke the quark-hadron duality prescription and replace 
the spectral density pf by 



Pi{p%s) 
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where ImFj (p^, s) is obtained from the imaginary part of the correlator functions calcu- 
lated in QCD starting from some threshold Sq. In order to suppress the higher states and 
continuum contributions, we follow the standard procedure in the QCD sum rules method 
and apply Borel transformation B on the variable (p + q)'^ to the dispersion integral to get, 



F,ip',M') = BF,{p\{p + qY) 

' p,{p\s)e~'"'"ds. (15) 



Using eqs.(10)-(12) and (14) we get 

F,{p\ M') = 2f'^{p')'^^^e'<"'" + - r ImF{3^^(p^ s)e-'/^"ds , (16) 

nib 71" "'so 
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+ - r ImF2«^^(p2, s)e-^/^"ds , (17) 

TT J so 

F,{p\M^) = _M^!<^e-™B/^' + - ^ImF3«^^(p^.)e-/^'rfs. (18) 
Mb + rriK rrib vr J so 

The main problem is then the calculation of the correlator functions (5) and (6) in QCD. 
After applying the Borel transformation, the result can be written in the following form: 

F,{p\ M') = - r Imi^^^^(/, s)e-'/^"ds . (19) 

TT Jml 



Equating the Fi in eq.(19) to the corresponding Fi in eqs.(16)-(18), we arrive at the sum 
rules for the formfactors, which describe the B ^ K transition: 



"b 
rso 



'2.'KfBm% 

rip') + rip') = ^^^ImF,^^^(p^.)e-(— i^)/^^d., (21) 

Mp') = ^^(mB + m;,) ^ImF3«^^(p^s)e-(— ^«)/^^rfs. (22) 
'^Jb'itT'B -'"T'b 

One can calculate ImFf^ (p"^, s), in the deep Euclidean region, where both p^ and (p + q)'^ 
are negative and large. The leading contribution to the operator product expansion comes 
from the contraction of the 6-quark operators expressed in eqs.(5) and (6) to the free 
6-quark propagator 

(o|rw.Wo)}|o>^->/|^e-i±^, 
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Note that, as can be seen from eqs.(23) and (24), the problem is reduced to the calculation 
of the matrix elements of the gauge-invariant, nonlocal operators sandwiched in between 
the vacuum and the K meson states. These matrix elements define the K meson light cone 
wave functions. Following [16-17], we define the i^-meson wave functions as: 



(ir(g)|s(x)7^(l-75)g(0)|0) = tq^K due 
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In eq.(25) (Pk{u) is the leading twist two, gi{u) and (72 (m) are the twist four K meson wave 
functions, respectively. The second matrix element of eq.(23), can be split into two matrix 
elements using the identity 7^7jy = g^u — '^o'^z/, and the result can be evaluated using the 
twist three wave functions defined as [16-19]: 
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The matrix elements in eq.(24) can be easily calculated using the identities. 
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to express it in terms of the wave functions in eq.(25). 

Substituting the matrix elements (25-27) into the eqs.(23) and (24) and integrating over 
the variables x and k we get the following expressions for Fi{p'^ , {p + q)'^)'- 
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where 
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, G2{u) = - g2{v)dv 
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There are also contributions to the above considered wave functions from muhi-particle 
meson wave functions. Here we consider only the operator qGq, which gives the main 
contribution and corresponds to the quark-antiquark-gluon components in the kaon (for 
more detail see pOf). In this approximation, the 6-quark propagator is defined as (see ||T8[| ): 
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It is clear that when we take into account the second term in the right hand side of eq.(31), 
new matrix elements appear. Substituting eq.(31) into eqs.(5) and (6), and using the 
identity 

IfMluC^pX = {(^fMxgup - (^ppQux) + i{9p\9up - 9i,p9u\) - 
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one can express the resulting new matrix elements in terms of the three particle wave 
functions [16-20]: 

{K{q)\s{x)gsG ^y{ux)ao,i3l^q{Q)\Gi) = if3K[iqf,qa9,yfS - qvqo,9t,p) - {.qtJiq(i9ua - quqi39fia) ] X 
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where Gap = ^taispxG''''^ , Vai = daida2das6{l — ai — a2 — as), and u = ai + uas. Here 
fsxio^i) is a twist three wave function and the remaining functions ip±, ip\\, (p± and (^y are 
all twist four wave functions. 

When we substitute eqs.(32)-(34) into eqs.(5) and (6), perform integration over x and 
k, add (28), (29), (30), apply the Borel transformation over {p + q)^ and equate the obtained 
results to (16), (17) and (18), we get the following sum rules for the B ^ K transition 
formfactors: 
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. The functions \E'±(\E'±), \1'||(\^||) in eqs.(35)-(37) are defined in the 
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Note that the formfactor f'^{p'^) is investigated in |]TB[ and [|r^ (In [|r^, it is necessary to 
make the simple replacement tt —>■ K), in light cone sum rules. Our results coincide with 
theirs if we let g^ = 0. 

At the end of this section we calculate the differential decay rate with the longitudinal 
polarization of the final leptons and we obtain 

+ 
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where A = 1 + r^ + s^ — 2r — 2s — 2sr, r = m\/m\, s = p^ /m\, ^ is the longitudinal 

polarization of the final lepton, m^ and v = xj "^ 

In eq.(38) A, C and D are defined as follows 



1 — —T- are its mass and velocity, respectively. 
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C = 2Cio/ ) 

D = Cio(/+ + /-) . 

For the dileptonic decays of the B mesons, the longitudinal polarization asymmetry Pl of 
the final state i, is defined by 

— (^ = -1)- — 
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(40) 



where ^ = — 1(+1) corresponds to the left (right) handed lepton. If in eq.(40), we let 
mi = 0, our results coincide with the results in [^ and if m£ ^ 0, they coincide with the 
ones in ""^^ 
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3 Numerical analysis 



The main input parameters in the sum rules (35)-(37) are the kaon wave functions on the 



light cone. For kaon wave functions we use the results of ref. |jT6|, |T^, ^ 

\2 
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6u(l — u) , 

2 ^ 
3 ^ 



u 



u){2u 



3 ' 



(p3K{ai) ~ 360a;ia2a: 
ip±{ai) ~ 105^(q;i - a2)«3 , 
(p\\{ai) ~ 1205^e{ai — a2)ttia2a3 , 
(f±{ai) ~ 105^03(1 - ag) , 
— 40(5^aia2«3 , 



m^^) 



here 6"^ (fib) — 0.17 GeV^ at fib — \/m% ~ '^l = 2.4 GeV"^, which follows from the QCD 
sum rules analysis (for more detail see |]18|), e{jib) — 0.36. In order to estimate hk, we use 
the PCAC relation for the pseudo Goldstone bosons 



/^x {{m) + {ss)) Pk 



fJ'TT 



^{qq)P 



0.62 (g = u or d) 



Here we use /^ ~ 133 MeV, fx — 160 MeV and we assume {qq) -^ {ss) = 1 ^ 0.8, which 
follows from QCD sum rules for strange hadrons [E^. As a result we have 



fiK — 1 GeV when /i^r 



m: 



uiu + md 



1.6 GeV . 



For the values of the other parameters, we choose: Jb — 0.14 GeV, which is obtained 
from 2-point QCD sum rules analysis [1^, ^, rUb — 4.7 GeV and Sq ^ 35 GeV^, and 
\VtbV:,\c:^QMb. 

Before giving numerical results on the formfactors, we must first determine the region 
for the Borel mass parameter M^, for which the sum rules yields reliable results. The 
lower limit of this region is determined by the requirement that, the terms proportional to 
M~'^"'{n > 1) remain subdominant. The upper limit of M^ is determined by requiring the 
higher resonance and continuum contributions to be less than ~ 30% of the total result. 
Our numerical results show that both requirements are satisfied in the region 8 GeV^ < 
M^ < 16 GeV^ and for the numerical analysis we use M^ = 10 GeV^. When p^ approaches 



the region ml — 0(1 GeV"^) a breakdown of the stabihty is expected, similar to the B ^ ir 
case (see |T|, [19| 



In fig.l, we present the p^ dependence of the formfactors f^ip^), f {p^) and frip^)- At 
zero momentum transfer, the QCD prediction for the formfactors are 

/+(p2 = 0) = 0.29 , 
/-(p2 = 0) = -0.21 , 
/t(/ = 0) = -0.31 . 

In fig.2a(b) we present the p^ dependence of the differential Branching ratios for the 
B -^ Kfi'^^~ and {B — >■ Kt^t^) decay with and without long distance effects. In both 
cases summation over the final lepton polarization is performed. Note that p^ dependence 



of the differential branching ratio B — > Kt^t is analysed in |11] and 0] using the light 



front formalism and heavy meson chiral theory, respectively. In [jlTI both short and long 



distance contributions are considered while in p4| only the short distance contributions are 



taken into account and in both works the obtained spectrum is fully symmetric while in 
the present work the spectrum we obtain seems to be slightly asymmetric as a result of 
a highly asymmetric resonance-type behaviour due to the nonperturbative contributions. 
Performing the integration over p^ in eq.(38) and using the values of the life times tb^ = 
(1.56 ±0.06) X 10~^^ s, tb^ = (1.62 ±0.06) x 10^^^ s ||2^, for the branching ratios, including 
only the short distance contributions, we get: 

B{Bd^ K^fi^fi') = (3.1±0.9) X 10"^ , 

B{Bd-^K\+T-) = (1.7±0.4) X 10^^ 

B{Bu^ K+fi^fi') = (3.2±0.8) X 10"^ , 

B{Bu -^ K+T+T-) = (1.77 ± 0.40) x 10"^ , 

where theoretical and experimental errors have been added quadratically. 

In fig. 3 we display the lepton longitudinal polarization asymmetry Pl as a function of 
p^ for the B -^ K^^ ^~ and B -^ Kt^t~ decays, at mt = 176 GeV, with and without the 
long distance contributions. From this figure one can see that Pl vanishes at the threshold 
due to the kinematical factor v and the value of Pl for the B -^ K^^^" decay varies in 
the region (0 ~. — 0.7) and (0 -. — 0.1) for the B -^ Kt^t~ decay, when long distance effects 
are excluded. 



4 Conclusion 

In this work, we calculate the transition formfactors for the exclusive B -^ Ki'^i~{i = fi, r) 
decay in the framework of the light cone QCD sum rules, and investigate the longitudinal 
polarization asymmetries of the muon and tau in this decay. From a comparison of our re- 
sults with the traditional QCD sum rule predictions (see 0), we observe that the behaviour 
of the formfactors are similar and the value of /^ in both approaches coincides, while /t 
differs two times than that of |^ at p^ = 0. It is important to note that for a more refined 
analysis, it is necessary to take into account the SU(3) breaking effects: the differences 
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between f^, hk and /^r, /Xtt and the differences of pion and kaon wave functions. Tliese 
SU(3) breaking terms can lead to differences between 5 — > vr and B ^ K formfactors. 
But we expect that these effects can change the results about 15-20% and this lies at the 
accuracy level of the sum rules method. 

Few words about the possibility of the experimental observation of this decay are in 
order. Experimentally, to observe an asymmetry Pi of a decay with the branching ratio 



B at the na level, the required number of events is A^ = -^-^ (see 111). For example, to 

observe the r lepton polarization at the exclusive channel B —>■ Kt^t^ at the 3a level, 
one needs at least A^ ~ 5 x 10^ B B decays. Since in the future 5-factories, it is expected 
that ~ 10^ S-mesons would be created per year, it is possible to measure the longitudinal 
polarization asymmetry of the r lepton. 
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Figure Captions 

1. The p^ dependence of the formfactors f^{p^), f'{p^) and frip"^)- 

2. a) Invariant mass squared distribution of the lepton pair for the decay B -^ Kfi^fi~. 
b) The same as in a), but for the decay B -^ Kt^t~ . 

Here and in all of the following figures the solid line corresponds to the short distance 
contributions only and the dashed line to the sum of both short and long distance contri- 
butions. 

3. a) The longitudinal polarization asymmetry P^ for the B -^ i^/i+/i~ decay, 
b) The same as in a), but for the B — > i^r"''r~ decay. 
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